3.95 \(\int \frac{(a+c x^2)^{3/2} (d+e x+f x^2)}{(g+h x)^4} \, dx\)

Optimal. Leaf size=475 \[ \frac{c \tanh ^{-1}\left (\frac{a h-c g x}{\sqrt{a+c x^2} \sqrt{a h^2+c g^2}}\right ) \left (3 a^2 h^4 (4 f g-e h)+3 a c g h^2 \left (11 f g^2-h (4 e g-d h)\right )+2 c^2 g^3 \left (10 f g^2-h (4 e g-d h)\right )\right )}{2 h^6 \left (a h^2+c g^2\right )^{3/2}}-\frac{\left (a+c x^2\right )^{5/2} \left (d h^2-e g h+f g^2\right )}{3 h (g+h x)^3 \left (a h^2+c g^2\right )}-\frac{\left (a+c x^2\right )^{3/2} \left (-x \left (3 a f h^2+c \left (5 f g^2-2 h (e g-d h)\right )\right )-3 a h (3 f g-e h)+c g \left (-d h+4 e g-\frac{10 f g^2}{h}\right )\right )}{6 h^2 (g+h x)^2 \left (a h^2+c g^2\right )}-\frac{\sqrt{a+c x^2} \left (c h x \left (3 a h^2 (3 f g-e h)+c g \left (10 f g^2-h (4 e g-d h)\right )\right )+\left (a h^2+c g^2\right ) \left (3 a f h^2+2 c \left (10 f g^2-h (4 e g-d h)\right )\right )\right )}{2 h^5 (g+h x) \left (a h^2+c g^2\right )}+\frac{\sqrt{c} \tanh ^{-1}\left (\frac{\sqrt{c} x}{\sqrt{a+c x^2}}\right ) \left (3 a f h^2+2 c \left (10 f g^2-h (4 e g-d h)\right )\right )}{2 h^6} \]

[Out]

-(((c*g^2 + a*h^2)*(3*a*f*h^2 + 2*c*(10*f*g^2 - h*(4*e*g - d*h))) + c*h*(3*a*h^2*(3*f*g - e*h) + c*g*(10*f*g^2
- h*(4*e*g - d*h)))*x)*Sqrt[a + c*x^2])/(2*h^5*(c*g^2 + a*h^2)*(g + h*x)) - ((c*g*(4*e*g - (10*f*g^2)/h - d*h
) - 3*a*h*(3*f*g - e*h) - (3*a*f*h^2 + c*(5*f*g^2 - 2*h*(e*g - d*h)))*x)*(a + c*x^2)^(3/2))/(6*h^2*(c*g^2 + a*
h^2)*(g + h*x)^2) - ((f*g^2 - e*g*h + d*h^2)*(a + c*x^2)^(5/2))/(3*h*(c*g^2 + a*h^2)*(g + h*x)^3) + (Sqrt[c]*(
3*a*f*h^2 + 2*c*(10*f*g^2 - h*(4*e*g - d*h)))*ArcTanh[(Sqrt[c]*x)/Sqrt[a + c*x^2]])/(2*h^6) + (c*(3*a^2*h^4*(4
*f*g - e*h) + 2*c^2*g^3*(10*f*g^2 - h*(4*e*g - d*h)) + 3*a*c*g*h^2*(11*f*g^2 - h*(4*e*g - d*h)))*ArcTanh[(a*h
- c*g*x)/(Sqrt[c*g^2 + a*h^2]*Sqrt[a + c*x^2])])/(2*h^6*(c*g^2 + a*h^2)^(3/2))

________________________________________________________________________________________

Rubi [A]  time = 0.844642, antiderivative size = 469, normalized size of antiderivative = 0.99, number of steps used = 8, number of rules used = 6, integrand size = 29, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.207, Rules used = {1651, 813, 844, 217, 206, 725} \[ \frac{c \tanh ^{-1}\left (\frac{a h-c g x}{\sqrt{a+c x^2} \sqrt{a h^2+c g^2}}\right ) \left (3 a^2 h^4 (4 f g-e h)+3 a c g h^2 \left (11 f g^2-h (4 e g-d h)\right )+2 c^2 \left (10 f g^5-g^3 h (4 e g-d h)\right )\right )}{2 h^6 \left (a h^2+c g^2\right )^{3/2}}-\frac{\left (a+c x^2\right )^{5/2} \left (d h^2-e g h+f g^2\right )}{3 h (g+h x)^3 \left (a h^2+c g^2\right )}-\frac{\left (a+c x^2\right )^{3/2} \left (-x \left (3 a f h^2-2 c h (e g-d h)+5 c f g^2\right )-3 a h (3 f g-e h)+c g \left (-d h+4 e g-\frac{10 f g^2}{h}\right )\right )}{6 h^2 (g+h x)^2 \left (a h^2+c g^2\right )}-\frac{\sqrt{a+c x^2} \left (c h x \left (3 a h^2 (3 f g-e h)-c g h (4 e g-d h)+10 c f g^3\right )+\left (a h^2+c g^2\right ) \left (3 a f h^2-2 c h (4 e g-d h)+20 c f g^2\right )\right )}{2 h^5 (g+h x) \left (a h^2+c g^2\right )}+\frac{\sqrt{c} \tanh ^{-1}\left (\frac{\sqrt{c} x}{\sqrt{a+c x^2}}\right ) \left (3 a f h^2-2 c h (4 e g-d h)+20 c f g^2\right )}{2 h^6} \]

Antiderivative was successfully verified.

[In]

Int[((a + c*x^2)^(3/2)*(d + e*x + f*x^2))/(g + h*x)^4,x]

[Out]

-(((c*g^2 + a*h^2)*(20*c*f*g^2 + 3*a*f*h^2 - 2*c*h*(4*e*g - d*h)) + c*h*(10*c*f*g^3 - c*g*h*(4*e*g - d*h) + 3*
a*h^2*(3*f*g - e*h))*x)*Sqrt[a + c*x^2])/(2*h^5*(c*g^2 + a*h^2)*(g + h*x)) - ((c*g*(4*e*g - (10*f*g^2)/h - d*h
) - 3*a*h*(3*f*g - e*h) - (5*c*f*g^2 + 3*a*f*h^2 - 2*c*h*(e*g - d*h))*x)*(a + c*x^2)^(3/2))/(6*h^2*(c*g^2 + a*
h^2)*(g + h*x)^2) - ((f*g^2 - e*g*h + d*h^2)*(a + c*x^2)^(5/2))/(3*h*(c*g^2 + a*h^2)*(g + h*x)^3) + (Sqrt[c]*(
20*c*f*g^2 + 3*a*f*h^2 - 2*c*h*(4*e*g - d*h))*ArcTanh[(Sqrt[c]*x)/Sqrt[a + c*x^2]])/(2*h^6) + (c*(3*a^2*h^4*(4
*f*g - e*h) + 3*a*c*g*h^2*(11*f*g^2 - h*(4*e*g - d*h)) + 2*c^2*(10*f*g^5 - g^3*h*(4*e*g - d*h)))*ArcTanh[(a*h
- c*g*x)/(Sqrt[c*g^2 + a*h^2]*Sqrt[a + c*x^2])])/(2*h^6*(c*g^2 + a*h^2)^(3/2))

Rule 1651

Int[(Pq_)*((d_) + (e_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> With[{Q = PolynomialQuotient[Pq, d
+ e*x, x], R = PolynomialRemainder[Pq, d + e*x, x]}, Simp[(e*R*(d + e*x)^(m + 1)*(a + c*x^2)^(p + 1))/((m + 1
)*(c*d^2 + a*e^2)), x] + Dist[1/((m + 1)*(c*d^2 + a*e^2)), Int[(d + e*x)^(m + 1)*(a + c*x^2)^p*ExpandToSum[(m
+ 1)*(c*d^2 + a*e^2)*Q + c*d*R*(m + 1) - c*e*R*(m + 2*p + 3)*x, x], x], x]] /; FreeQ[{a, c, d, e, p}, x] && Po
lyQ[Pq, x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[m, -1]

Rule 813

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Simp[((d + e*x)^(
m + 1)*(e*f*(m + 2*p + 2) - d*g*(2*p + 1) + e*g*(m + 1)*x)*(a + c*x^2)^p)/(e^2*(m + 1)*(m + 2*p + 2)), x] + Di
st[p/(e^2*(m + 1)*(m + 2*p + 2)), Int[(d + e*x)^(m + 1)*(a + c*x^2)^(p - 1)*Simp[g*(2*a*e + 2*a*e*m) + (g*(2*c
*d + 4*c*d*p) - 2*c*e*f*(m + 2*p + 2))*x, x], x], x] /; FreeQ[{a, c, d, e, f, g, m}, x] && NeQ[c*d^2 + a*e^2,
0] && RationalQ[p] && p > 0 && (LtQ[m, -1] || EqQ[p, 1] || (IntegerQ[p] &&  !RationalQ[m])) && NeQ[m, -1] &&
!ILtQ[m + 2*p + 1, 0] && (IntegerQ[m] || IntegerQ[p] || IntegersQ[2*m, 2*p])

Rule 844

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Dist[g/e, Int[(d
+ e*x)^(m + 1)*(a + c*x^2)^p, x], x] + Dist[(e*f - d*g)/e, Int[(d + e*x)^m*(a + c*x^2)^p, x], x] /; FreeQ[{a,
c, d, e, f, g, m, p}, x] && NeQ[c*d^2 + a*e^2, 0] &&  !IGtQ[m, 0]

Rule 217

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Subst[Int[1/(1 - b*x^2), x], x, x/Sqrt[a + b*x^2]] /; FreeQ[{a,
b}, x] &&  !GtQ[a, 0]

Rule 206

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTanh[(Rt[-b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[-b, 2]), x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 725

Int[1/(((d_) + (e_.)*(x_))*Sqrt[(a_) + (c_.)*(x_)^2]), x_Symbol] :> -Subst[Int[1/(c*d^2 + a*e^2 - x^2), x], x,
(a*e - c*d*x)/Sqrt[a + c*x^2]] /; FreeQ[{a, c, d, e}, x]

Rubi steps

\begin{align*} \int \frac{\left (a+c x^2\right )^{3/2} \left (d+e x+f x^2\right )}{(g+h x)^4} \, dx &=-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}-\frac{\int \frac{\left (-3 (c d g-a f g+a e h)-\left (3 a f h-c \left (2 e g-\frac{5 f g^2}{h}-2 d h\right )\right ) x\right ) \left (a+c x^2\right )^{3/2}}{(g+h x)^3} \, dx}{3 \left (c g^2+a h^2\right )}\\ &=-\frac{\left (c g \left (4 e g-\frac{10 f g^2}{h}-d h\right )-3 a h (3 f g-e h)-\left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right ) x\right ) \left (a+c x^2\right )^{3/2}}{6 h^2 \left (c g^2+a h^2\right ) (g+h x)^2}-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}+\frac{\int \frac{\left (4 a \left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right )-\frac{4 c \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right ) x}{h}\right ) \sqrt{a+c x^2}}{(g+h x)^2} \, dx}{8 h^2 \left (c g^2+a h^2\right )}\\ &=-\frac{\left (\left (c g^2+a h^2\right ) \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )+c h \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right ) x\right ) \sqrt{a+c x^2}}{2 h^5 \left (c g^2+a h^2\right ) (g+h x)}-\frac{\left (c g \left (4 e g-\frac{10 f g^2}{h}-d h\right )-3 a h (3 f g-e h)-\left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right ) x\right ) \left (a+c x^2\right )^{3/2}}{6 h^2 \left (c g^2+a h^2\right ) (g+h x)^2}-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}-\frac{\int \frac{8 a c \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right )-\frac{8 c \left (c g^2+a h^2\right ) \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right ) x}{h}}{(g+h x) \sqrt{a+c x^2}} \, dx}{16 h^4 \left (c g^2+a h^2\right )}\\ &=-\frac{\left (\left (c g^2+a h^2\right ) \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )+c h \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right ) x\right ) \sqrt{a+c x^2}}{2 h^5 \left (c g^2+a h^2\right ) (g+h x)}-\frac{\left (c g \left (4 e g-\frac{10 f g^2}{h}-d h\right )-3 a h (3 f g-e h)-\left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right ) x\right ) \left (a+c x^2\right )^{3/2}}{6 h^2 \left (c g^2+a h^2\right ) (g+h x)^2}-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}+\frac{\left (c \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )\right ) \int \frac{1}{\sqrt{a+c x^2}} \, dx}{2 h^6}-\frac{\left (c \left (3 a^2 h^4 (4 f g-e h)+3 a c g h^2 \left (11 f g^2-h (4 e g-d h)\right )+2 c^2 \left (10 f g^5-g^3 h (4 e g-d h)\right )\right )\right ) \int \frac{1}{(g+h x) \sqrt{a+c x^2}} \, dx}{2 h^6 \left (c g^2+a h^2\right )}\\ &=-\frac{\left (\left (c g^2+a h^2\right ) \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )+c h \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right ) x\right ) \sqrt{a+c x^2}}{2 h^5 \left (c g^2+a h^2\right ) (g+h x)}-\frac{\left (c g \left (4 e g-\frac{10 f g^2}{h}-d h\right )-3 a h (3 f g-e h)-\left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right ) x\right ) \left (a+c x^2\right )^{3/2}}{6 h^2 \left (c g^2+a h^2\right ) (g+h x)^2}-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}+\frac{\left (c \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )\right ) \operatorname{Subst}\left (\int \frac{1}{1-c x^2} \, dx,x,\frac{x}{\sqrt{a+c x^2}}\right )}{2 h^6}+\frac{\left (c \left (3 a^2 h^4 (4 f g-e h)+3 a c g h^2 \left (11 f g^2-h (4 e g-d h)\right )+2 c^2 \left (10 f g^5-g^3 h (4 e g-d h)\right )\right )\right ) \operatorname{Subst}\left (\int \frac{1}{c g^2+a h^2-x^2} \, dx,x,\frac{a h-c g x}{\sqrt{a+c x^2}}\right )}{2 h^6 \left (c g^2+a h^2\right )}\\ &=-\frac{\left (\left (c g^2+a h^2\right ) \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right )+c h \left (10 c f g^3-c g h (4 e g-d h)+3 a h^2 (3 f g-e h)\right ) x\right ) \sqrt{a+c x^2}}{2 h^5 \left (c g^2+a h^2\right ) (g+h x)}-\frac{\left (c g \left (4 e g-\frac{10 f g^2}{h}-d h\right )-3 a h (3 f g-e h)-\left (5 c f g^2+3 a f h^2-2 c h (e g-d h)\right ) x\right ) \left (a+c x^2\right )^{3/2}}{6 h^2 \left (c g^2+a h^2\right ) (g+h x)^2}-\frac{\left (f g^2-e g h+d h^2\right ) \left (a+c x^2\right )^{5/2}}{3 h \left (c g^2+a h^2\right ) (g+h x)^3}+\frac{\sqrt{c} \left (20 c f g^2+3 a f h^2-2 c h (4 e g-d h)\right ) \tanh ^{-1}\left (\frac{\sqrt{c} x}{\sqrt{a+c x^2}}\right )}{2 h^6}+\frac{c \left (3 a^2 h^4 (4 f g-e h)+3 a c g h^2 \left (11 f g^2-h (4 e g-d h)\right )+2 c^2 \left (10 f g^5-g^3 h (4 e g-d h)\right )\right ) \tanh ^{-1}\left (\frac{a h-c g x}{\sqrt{c g^2+a h^2} \sqrt{a+c x^2}}\right )}{2 h^6 \left (c g^2+a h^2\right )^{3/2}}\\ \end{align*}

Mathematica [A]  time = 1.40119, size = 517, normalized size = 1.09 \[ \frac{-\frac{h \sqrt{a+c x^2} \left ((g+h x)^2 \left (6 a^2 f h^4+a c h^2 \left (h (8 d h-23 e g)+50 f g^2\right )+c^2 \left (g^2 h (11 d h-26 e g)+47 f g^4\right )\right )-(g+h x) \left (a h^2+c g^2\right ) \left (-3 a h^2 (e h-2 f g)+c g h (7 d h-10 e g)+13 c f g^3\right )+2 \left (a h^2+c g^2\right )^2 \left (h (d h-e g)+f g^2\right )+6 c (g+h x)^3 \left (a h^2+c g^2\right ) (4 f g-e h)-3 c f h x (g+h x)^3 \left (a h^2+c g^2\right )\right )}{(g+h x)^3 \left (a h^2+c g^2\right )}+\frac{3 c \log \left (\sqrt{a+c x^2} \sqrt{a h^2+c g^2}+a h-c g x\right ) \left (-3 a^2 h^4 (e h-4 f g)+3 a c g h^2 \left (h (d h-4 e g)+11 f g^2\right )+2 c^2 \left (g^3 h (d h-4 e g)+10 f g^5\right )\right )}{\left (a h^2+c g^2\right )^{3/2}}-\frac{3 c \log (g+h x) \left (-3 a^2 h^4 (e h-4 f g)+3 a c g h^2 \left (h (d h-4 e g)+11 f g^2\right )+2 c^2 \left (g^3 h (d h-4 e g)+10 f g^5\right )\right )}{\left (a h^2+c g^2\right )^{3/2}}+3 \sqrt{c} \log \left (\sqrt{c} \sqrt{a+c x^2}+c x\right ) \left (3 a f h^2+2 c h (d h-4 e g)+20 c f g^2\right )}{6 h^6} \]

Antiderivative was successfully verified.

[In]

Integrate[((a + c*x^2)^(3/2)*(d + e*x + f*x^2))/(g + h*x)^4,x]

[Out]

(-((h*Sqrt[a + c*x^2]*(2*(c*g^2 + a*h^2)^2*(f*g^2 + h*(-(e*g) + d*h)) - (c*g^2 + a*h^2)*(13*c*f*g^3 + c*g*h*(-
10*e*g + 7*d*h) - 3*a*h^2*(-2*f*g + e*h))*(g + h*x) + (6*a^2*f*h^4 + a*c*h^2*(50*f*g^2 + h*(-23*e*g + 8*d*h))
+ c^2*(47*f*g^4 + g^2*h*(-26*e*g + 11*d*h)))*(g + h*x)^2 + 6*c*(4*f*g - e*h)*(c*g^2 + a*h^2)*(g + h*x)^3 - 3*c
*f*h*(c*g^2 + a*h^2)*x*(g + h*x)^3))/((c*g^2 + a*h^2)*(g + h*x)^3)) - (3*c*(-3*a^2*h^4*(-4*f*g + e*h) + 3*a*c*
g*h^2*(11*f*g^2 + h*(-4*e*g + d*h)) + 2*c^2*(10*f*g^5 + g^3*h*(-4*e*g + d*h)))*Log[g + h*x])/(c*g^2 + a*h^2)^(
3/2) + 3*Sqrt[c]*(20*c*f*g^2 + 3*a*f*h^2 + 2*c*h*(-4*e*g + d*h))*Log[c*x + Sqrt[c]*Sqrt[a + c*x^2]] + (3*c*(-3
*a^2*h^4*(-4*f*g + e*h) + 3*a*c*g*h^2*(11*f*g^2 + h*(-4*e*g + d*h)) + 2*c^2*(10*f*g^5 + g^3*h*(-4*e*g + d*h)))
*Log[a*h - c*g*x + Sqrt[c*g^2 + a*h^2]*Sqrt[a + c*x^2]])/(c*g^2 + a*h^2)^(3/2))/(6*h^6)

________________________________________________________________________________________

Maple [B]  time = 0.252, size = 9835, normalized size = 20.7 \begin{align*} \text{output too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((c*x^2+a)^(3/2)*(f*x^2+e*x+d)/(h*x+g)^4,x)

[Out]

result too large to display

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c*x^2+a)^(3/2)*(f*x^2+e*x+d)/(h*x+g)^4,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c*x^2+a)^(3/2)*(f*x^2+e*x+d)/(h*x+g)^4,x, algorithm="fricas")

[Out]

Timed out

________________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{\left (a + c x^{2}\right )^{\frac{3}{2}} \left (d + e x + f x^{2}\right )}{\left (g + h x\right )^{4}}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c*x**2+a)**(3/2)*(f*x**2+e*x+d)/(h*x+g)**4,x)

[Out]

Integral((a + c*x**2)**(3/2)*(d + e*x + f*x**2)/(g + h*x)**4, x)

________________________________________________________________________________________

Giac [B]  time = 1.68426, size = 2565, normalized size = 5.4 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c*x^2+a)^(3/2)*(f*x^2+e*x+d)/(h*x+g)^4,x, algorithm="giac")

[Out]

1/2*sqrt(c*x^2 + a)*(c*f*x/h^4 - 2*(4*c*f*g*h^10 - c*h^11*e)/h^15) - (20*c^3*f*g^5 + 2*c^3*d*g^3*h^2 + 33*a*c^
2*f*g^3*h^2 + 3*a*c^2*d*g*h^4 + 12*a^2*c*f*g*h^4 - 8*c^3*g^4*h*e - 12*a*c^2*g^2*h^3*e - 3*a^2*c*h^5*e)*arctan(
-((sqrt(c)*x - sqrt(c*x^2 + a))*h + sqrt(c)*g)/sqrt(-c*g^2 - a*h^2))/((c*g^2*h^6 + a*h^8)*sqrt(-c*g^2 - a*h^2)
) - 1/3*(60*(sqrt(c)*x - sqrt(c*x^2 + a))^5*c^3*f*g^5*h^2 + 18*(sqrt(c)*x - sqrt(c*x^2 + a))^5*c^3*d*g^3*h^4 +
69*(sqrt(c)*x - sqrt(c*x^2 + a))^5*a*c^2*f*g^3*h^4 + 15*(sqrt(c)*x - sqrt(c*x^2 + a))^5*a*c^2*d*g*h^6 + 12*(s
qrt(c)*x - sqrt(c*x^2 + a))^5*a^2*c*f*g*h^6 - 36*(sqrt(c)*x - sqrt(c*x^2 + a))^5*c^3*g^4*h^3*e - 36*(sqrt(c)*x
- sqrt(c*x^2 + a))^5*a*c^2*g^2*h^5*e - 3*(sqrt(c)*x - sqrt(c*x^2 + a))^5*a^2*c*h^7*e + 210*(sqrt(c)*x - sqrt(
c*x^2 + a))^4*c^(7/2)*f*g^6*h + 54*(sqrt(c)*x - sqrt(c*x^2 + a))^4*c^(7/2)*d*g^4*h^3 + 183*(sqrt(c)*x - sqrt(c
*x^2 + a))^4*a*c^(5/2)*f*g^4*h^3 + 27*(sqrt(c)*x - sqrt(c*x^2 + a))^4*a*c^(5/2)*d*g^2*h^5 - 18*(sqrt(c)*x - sq
rt(c*x^2 + a))^4*a^2*c^(3/2)*f*g^2*h^5 - 12*(sqrt(c)*x - sqrt(c*x^2 + a))^4*a^2*c^(3/2)*d*h^7 - 6*(sqrt(c)*x -
sqrt(c*x^2 + a))^4*a^3*sqrt(c)*f*h^7 - 120*(sqrt(c)*x - sqrt(c*x^2 + a))^4*c^(7/2)*g^5*h^2*e - 84*(sqrt(c)*x
- sqrt(c*x^2 + a))^4*a*c^(5/2)*g^3*h^4*e + 21*(sqrt(c)*x - sqrt(c*x^2 + a))^4*a^2*c^(3/2)*g*h^6*e + 188*(sqrt(
c)*x - sqrt(c*x^2 + a))^3*c^4*f*g^7 + 44*(sqrt(c)*x - sqrt(c*x^2 + a))^3*c^4*d*g^5*h^2 - 82*(sqrt(c)*x - sqrt(
c*x^2 + a))^3*a*c^3*f*g^5*h^2 - 34*(sqrt(c)*x - sqrt(c*x^2 + a))^3*a*c^3*d*g^3*h^4 - 276*(sqrt(c)*x - sqrt(c*x
^2 + a))^3*a^2*c^2*f*g^3*h^4 - 48*(sqrt(c)*x - sqrt(c*x^2 + a))^3*a^2*c^2*d*g*h^6 - 36*(sqrt(c)*x - sqrt(c*x^2
+ a))^3*a^3*c*f*g*h^6 - 104*(sqrt(c)*x - sqrt(c*x^2 + a))^3*c^4*g^6*h*e + 64*(sqrt(c)*x - sqrt(c*x^2 + a))^3*
a*c^3*g^4*h^3*e + 138*(sqrt(c)*x - sqrt(c*x^2 + a))^3*a^2*c^2*g^2*h^5*e - 354*(sqrt(c)*x - sqrt(c*x^2 + a))^2*
a*c^(7/2)*f*g^6*h - 78*(sqrt(c)*x - sqrt(c*x^2 + a))^2*a*c^(7/2)*d*g^4*h^3 - 276*(sqrt(c)*x - sqrt(c*x^2 + a))
^2*a^2*c^(5/2)*f*g^4*h^3 - 36*(sqrt(c)*x - sqrt(c*x^2 + a))^2*a^2*c^(5/2)*d*g^2*h^5 + 60*(sqrt(c)*x - sqrt(c*x
^2 + a))^2*a^3*c^(3/2)*f*g^2*h^5 + 12*(sqrt(c)*x - sqrt(c*x^2 + a))^2*a^3*c^(3/2)*d*h^7 + 12*(sqrt(c)*x - sqrt
(c*x^2 + a))^2*a^4*sqrt(c)*f*h^7 + 192*(sqrt(c)*x - sqrt(c*x^2 + a))^2*a*c^(7/2)*g^5*h^2*e + 114*(sqrt(c)*x -
sqrt(c*x^2 + a))^2*a^2*c^(5/2)*g^3*h^4*e - 48*(sqrt(c)*x - sqrt(c*x^2 + a))^2*a^3*c^(3/2)*g*h^6*e + 222*(sqrt(
c)*x - sqrt(c*x^2 + a))*a^2*c^3*f*g^5*h^2 + 48*(sqrt(c)*x - sqrt(c*x^2 + a))*a^2*c^3*d*g^3*h^4 + 231*(sqrt(c)*
x - sqrt(c*x^2 + a))*a^3*c^2*f*g^3*h^4 + 33*(sqrt(c)*x - sqrt(c*x^2 + a))*a^3*c^2*d*g*h^6 + 24*(sqrt(c)*x - sq
rt(c*x^2 + a))*a^4*c*f*g*h^6 - 120*(sqrt(c)*x - sqrt(c*x^2 + a))*a^2*c^3*g^4*h^3*e - 102*(sqrt(c)*x - sqrt(c*x
^2 + a))*a^3*c^2*g^2*h^5*e + 3*(sqrt(c)*x - sqrt(c*x^2 + a))*a^4*c*h^7*e - 47*a^3*c^(5/2)*f*g^4*h^3 - 11*a^3*c
^(5/2)*d*g^2*h^5 - 50*a^4*c^(3/2)*f*g^2*h^5 - 8*a^4*c^(3/2)*d*h^7 - 6*a^5*sqrt(c)*f*h^7 + 26*a^3*c^(5/2)*g^3*h
^4*e + 23*a^4*c^(3/2)*g*h^6*e)/((c*g^2*h^6 + a*h^8)*((sqrt(c)*x - sqrt(c*x^2 + a))^2*h + 2*(sqrt(c)*x - sqrt(c
*x^2 + a))*sqrt(c)*g - a*h)^3) - 1/2*(20*c^(3/2)*f*g^2 + 2*c^(3/2)*d*h^2 + 3*a*sqrt(c)*f*h^2 - 8*c^(3/2)*g*h*e
)*log(abs(-sqrt(c)*x + sqrt(c*x^2 + a)))/h^6